AD-A231  327 


OFFICE  OF  NAVAL  RESEARCH 


OTIC  riLE  COPY 


Grant  N00014-90-J-1193 
TECHNICAL  REPORT  No.  41 

Critical  Behavior  in  Magnetic  Superlattices 
by 

T.  Hai,  Z.  Y.  Li,  D.  L.  Lin  and  Thomas  F.  George 


Prepared  for  publication 
in 

Journal  of  Magnetism  and  Magnetic  Materials 

Departments  of  Chemistry  and  Physics 
State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


DTIC 

I  electe 
JAN  3  01991 


January  1991 


Reproduction  in  whole  or  in  pan  is  permitted  for  any  purpose  of  the 
United  States  Government. 

This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /  DOWNGRADING  SCHCDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 
UBUFFALO/DC/9 1 /TR-4 1 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approval 
OMB  NO.  0704-0188 


6b.  OFFICE  SYMBOL 
(If  applicable) 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 

unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


8b.  OFFICE  SYMBOL 
(If  applicant) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Depts.  Chemistry  &  Physics 
State  University  of  New  York 


6c  AOORESS  (City,  State,  and  ZIP  Coda) 

Fronczak  Hall,  Amherst  Campus 
Buffalo,  New  York  14260 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

Office  of  Naval  Research  -  . 


8c  ADDRESS  (City,  State,  and  ZIP  Coda) 
Chemistry  Program 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


1 1 .  TITLE  (Include  Security  Classification) 

Critical  Behavior  in  Magnetic  Superlattices 


12.  PERSONAL  AUTHORS) 

T.  Hai,  Z.  Y.  Li,  D.  L.  Lin  and  Thomas  F.  George 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Year,  Month,  Day)  IS.  PAGE  COUNT 

from _ to _  January  1991 _ 20 


16.  SUPPLEMENTARY  NOTATION 

Prepared  for  publication 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Chemistry  Program 

800  N.  Quincy  Street 

Arlington,  Virginia  22217  _ 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Grant  N00014-90-J-1193 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

element  no. 


PROJECT 

TASK 

NO. 

NO 

Journal  of  Magnetism  and  Magnetic  Materials 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

earnGAL  hhiaviCr*  gorrclatioho  - 

-MAGNETIC  SUPERLATTICES,  FERROMAGNETS  j£r~ 

^EFFECTIVE  FIELD  THEORY  CIIMfi"«ggfg  Mgggfr _ 


RACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


Based  on  the  effective  field  theory  with  correlations,  the  critical  behavior  of  a 
magnetic  superlattice  consisting  of  two  different  ferromagnets  is  examined.  A  simple 
cubic  Ising  model  with  nearest-neighbor  coupling  is  assumed.  It  is  found  that  there 
exists  a  critical  value  for  the  interface  exchange  coupling  above  which  the  interface 
magnetism  appears.  The  critical  temperature  of  the  superlattice  is  studied  as  a  function 
of  the  thickness  of  Che  constituents  in  a  unit  cell  and  of  the  exchange  coupling  energies 
in  each  material.  A  phase  diagram  is  also  given. 

4 


20.  DISTRIBUTION  /AVAILABILITY  OF  ABSTRACT 

Bunclassifieo/unumited  B  same  as  rpt 


224.  NAME  OF  RESPONSIBLE  INDIVIDUAL 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

dtic  users  I  Unclassified 

22b  TELEPHONE  (Include  Area  Code)  I  22c.  office  symbol 


r<lr4M-vI.T3T" 


DO  Form  1473,  JUN  M 


Previous  editions  are  obsolete. 


UNCLASSIFIED 


Journal  of  Magnetism  and  Magnetic  Materials,  in  press 


Critical  behavior  in  magnetic  superlattices 

T.  Hai 

P.  0.  Box  1001-73 
Zhengzhou  450002,  P.  R.  China 

* 

Z.  Y.  Li  ,  D.  L.  Lin  and  Thomas  F.  George 
Department  of  Physics  and  Astronomy 
State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260,  USA 

Abstract 

Based  on  the  effective  field  theory  with  correlations,  the  critical 
behavior  of  a  magnetic  superlattice  consisting  of  two  different  ferromagnets 
is  examined.  A  simple  cubic  Ising  model  with  nearest -neighbor  coupling  is 
assumed.  It  is  found  that  there  exists  a  critical  value  for  the  interface 
exchange  coupling  above  which  the  interface  magnetism  appears.  The  critical 
temperature  of  the  superlattice  is  studied  as  a  function  of  the  thickness  of 
the  constituents  in  a  unit  cell  and  of  the  exchange  coupling  energies  in  each 
material.  A  phase  diagram  is  also  given. 


*  On  leave  of  absence  from  the  Department  of  Physics,  Suzhou  University, 
Suzhou  215006,  P.  R.  China 


Acce;io.i  For  ^ 

NTiS 

CfM&l  ^ 

DTiC 

'1 A  J  □ 

U:  a.’i : 

oii  :cj  |_i 

JilUliiCUtiOii 

%7 

Diit  ibution  / 

Availability  Codes 

Dist 

Avail  and  lor 

Special 

2 


I.  Introduction 

In  recent  years,  there  has  been  increasing  interest  in  the  nature  of 

spin  waves  as  well  as  critical  phenomena  in  various  magnetic  layered 

structures  and  superlattices.  Ha  and  Tsai^  have  studied  the  variation  with 

modulation  wavelength  of  the  Curie  temperature  for  a  Heisenberg  magnetic 

superlattice.  Their  results  agree  qualitatively  with  experiments  on  the  Cu/Ni 
2 

film.  Superlattice  structures  composed  of  alternating  ferromagnetic  and 

3  4 

antiferromagnetic  layers  have  been  investigated  by  Hinchey  and  Mills,  ’  using 
a  localized  spin  model.  A  sequence  of  spin- reorientation  transitions  are 
found  to  be  different  for  superlattices  with  the  antiferromagnetic  component 
consisting  of  an  even  or  odd  number  of  spin  layers.  In  addition,  the  spin- 
wave  spectrum  and  infrared  absorption  spectrum  are  calculated. 

For  a  periodic  multilayer  system  formed  from  two  different  ferromagnetic 
materials,  Fishman  et  al5  have  discussed  its  statics  and  dynamics  within  the 
framework  of  the  Ginzburg- Landau  formulation.  They  have  computed  the 
transition  temperature  and  spin-wave  spectra.  On  the  other  hand,  the  Landau 
formalism  of  Camley  and  Tilley^  has  been  applied  to  calculate  the  critical 
temperature  in  the  same  system.^  Compared  to  Ref.  5,  the  formalism  of  Ref.  6 
appears  to  be  more  general  because  it  allows  for  a  wider  range  of  boundary 
conditions  and  includes  the  sign  of  exchange  coupling  across  the  interface. 

For  more  complicated  superlattices  with  arbitrary  number  of  different 

g 

layers  in  an  elementary  unit,  Barnas  has  derived  some  general  dispersion 
equations  for  the  bulk  and  surface  magnetic  polaritons.  These  equations  are 
then  applied  to  magnetostatic  modes  and  to  retarded  wave  propagation  in  the 

9 

Voigt  geometry. 

We  study,  in  this  article,  the  critical  temperature  in  an  infinite 
superlattice  consisting  of  two  ferromagnetic  materials  with  different  bulk 


3 


properties.  In  particular,  we  consider  the  two  constituents  A  and  B  with 

different  bulk  transition  temperatures  as  a  simple  model,  i.e.,  T^a^  *  . 

c  c 

The  interface  is  in  general  different  in  nature  from  both  bulks,  even  if  the 
bulk  critical  temperatures  are  the  same.  We  use  the  effective -field  theory 
with  correlations^  in  the  present  work,  as  it  is  believed  to  be  far  superior 
to  the  standard  mean- field  approximation. 

Because  of  the  periodicity  of  the  superlattice  structure,  we  restrict 
our  discussions  to  a  unit  cell  which  interacts  with  its  nearest-neighbor  cells 
via  interface  couplings.  Thus  the  system  can  be  treated  by  extending  the 
method  developed  for  a  magnetic  slab.^  Our  major  concerns  are  the  dependence 
of  the  transition  temperature  on  the  thickness  of  individual  constituents  in 
the  cell  and  the  influence  of  the  interface  magnetic  properties  on  the  phase 
transition  temperature.  These  questions,  to  our  knowledge,  have  not  been 
considered  in  the  literature  thus  far. 

In  Sec.  II.  we  outline  the  theory  and  derive  the  equation  that 
determines  the  transition  temperature.  Numerical  results  are  discussed  in 
Sec.  Ill  where  the  existence  of  the  interface  magnetic  phase  transition  is 
discovered  and  the  critical  value  of  the  interface  coupling  relative  the  bulk 
coupling  is  determined.  A  brief  conclusion  is  given  in  Sec.  IV. 

II .  Theory 

Consider  an  infinite  superlattice  consisting  of  two  different 
ferromagnetic  materials  A  and  B.  For  simplicity,  we  restrict  our  attention  to 
the  case  of  simple  cubic  Ising-type  structures.  The  periodic  condition 
suggests  that  we  only  have  to  consider  one  unit  cell  which  interacts  with  its 
nearest  neighbors  via  the  interface  coupling.  The  situation  is  depicted  in 
Fig.  1.  The  coupling  strength  between  nearest-neighboring  spins  in  A(B)  is 


4 


denoted  by  ,  while  stands  for  the  exchange  coupling  between  the 
nearest-neighbor  spins  across  the  interface.  The  corresponding  number  of 
atomic  layers  in  A(B)  is  1^(1^) ,  and  the  thickness  of  the  cell  is  L  -  Lfl  +  L^. 
The  Hamiltonian  of  the  system  is  given  by 


H-  ’  2^ 

l.j 


Visj 


where  the  sum  is  taken  over  all  the  nearest-neighbor  pairs  only  once,  -  ±1 

is  the  usual  Ising  variable,  and  stands  for  one  of  the  three  coupling 

constants  depending  on  where  the  spin  pair  is  located. 

To  evaluate  the  mean  spin  <S^>,  we  start  with  the  exact  Callen 
12 

identity, 


<Sj>  -  <tanh  (/3 


}  Jii  sj 


)>  . 


where  /?  -  l/k_T,  <...>  indicates  the  usual  canonical  ensemble  average  for  a 

D 

given  configuration  of  { J ^ ^ ) ,  and  j  runs  over  all  nearest  neighbors  of  site  i. 

We  now  introduce  the  differential  operator  D  -  and  recall  that  the 

ox 

displacement  operator  exp(aO)  satisfies  the  relation 


exp(aD)f(x)  -  f(x+a) 


Equation  (2)  can  then  be  rewritten,  with  the  help  of  (3),  as 


<Sj>  -  <exp(D  ^  [tanh(/Jx)l: 

j 
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<|ll  [cosh(DJ.^)  +  Sj  sinh(DJij)]>  tanh(^x)  |x_q 


This  equaCion  represents  a  set  of  L  coupled  equations  for  the  L  spin  layers  in 

the  unit  cell.  Each  layer  is  only  coupled  to  its  nearest-neighbor  layers. 

The  multi-spin  correlation  function,  however,  must  be  decoupled  before  a 

13 

practical  calculation  can  be  made.  We  follow  the  standard  procedure  with 
the  decoupling  approximation  <x^x2..x^>  -  <x^Xx^>  ...  <x^>.  Thus  the  average 
atomic  magnetization  of  each  layer  is  given  by 


-  [cosh(DJa^)  +  mosinh(DJab) ] [cosh(DJa>  +  m^sinhCDJ^) 


x  [cosh(DJ  )  +  ra,sinh(DJ  )]f(x)  n 

3  c.  3  X— U 


m^  -  (cosh(DJa)  +  m.  ^sinh(DJa> } (cosh(DJa)  +  nusinh(DJa)] 


x  (cosh(DJa)  +  m.+1sinh(DJa> ] f (x)  2  <  i  <  Lg 


dl  -  (cosh(DJ  )  +  m  ^inMOJ  )][cosh(DJ  )  +  mT  sinh(DJa)] 
L  3  L  -  I  3  3  L 


x  (cosh(DJab)  +  ^  +Isinh(DJab)]f(x)  x_Q 


-  [cosh(DJflb)  +  m^  sinh(DJab> ] [cosh(DJb>  +  +^sinh(DJb>] 


x  [cosh(DJb)  +  +2sinh(DJab)}f(x)  x_0 
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4 

m.  -  [cosh(DJ  )  +  m,  . sinh(DJ,  ) ] [cosh(DJ.  )  +  m,sinh(DJ,  )] 

J  aj-1  b  b  j  b 

[cosh(DJb)  +  mj+1sinh(DJb)]f(x)  x-Q  <  j  <  L-l  (5e) 

4 

-  [cosh(DJb>  +  ^sinh(DJb) ]cosh(DJb)  +  mr sinh(DJb> ] 

[cosh(1XJab)  +  mL+lsinh<DJab)lf(x)|x-0  ’  (5f) 

where  we  have  defined 

f(x)  -  tanh(^x)  (6) 

on  which  the  operator  D  applies.  The  variable  x  is  set  to  zero  at  the  end  of 
the  calculation  in  each  equation.  It  is  observed  that  the  periodic  condition 
of  the  superlattice  is  indeed  satisfied,  namely,  mQ  -  m^  and  -  m^^.  As 
the  temperature  becomes  higher  than  the  critical  temperature,  the  whole  system 
becomes  demagnetized.  Thus,  we  can  determine  Tc  from  Eq.  (5)  by  requiring 
that  the  mean  atomic  magnetization  in  every  spin  layer  approaches  zero. 
Consequently,  all  terms  of  the  order  higher  than  linear  in  Eqs.  (5)  can  be 
neglected.  This  leads  to  a  set  of  simultaneous  equations 

(l-4A1)m1  -  A2mQ  -  A^s,  -  0 

<l-4A3)m1  -  A3 (“i-l+mi+i>  “  0 
# 

U-4A1>\.  •  Vl  -1  •  *2*1.  *1  * 

a  a  a 


0 


0 


(l-4A4)m^  +1  *  A5mL  +2 

a  a  a 

<l-4A6)“j  *  A6<mj.l-mj+l)  “  ° 

<l-4A4>n,L  *  Vl-1  ‘  A5mL+l  “  °  * 

where  the  coefficients  are  given  by 

4 

-  cosh  (DJa)sinh(DJa)cosh(DJab)f (x)  x_q 

A2  -  cosh5(DJa)sinh(DJab)f(x)  x_Q 

A,  -  cosh^ (DJ  )sinh(DJ  )f(x)  n 
J  a  a  x**U 

A^  -  cosh^(DJb)sinh(DJb)cosh(DJab)f(x)  x_q 
A5  -  cosh5(DJb)sinh(DJab)f(x)  x-Q 
A6  -  cosh5(DJb)sinh(DJb)f<x)  x-Q 


The  secular  equation  of  this  set  of  coupled  equations  is 
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where  we  have  defined 

ax  -  (1-4A1)/A1  a2  -  (1-4A4)/A4 

bl  "  ‘A2/Al  b2  “  'A5/A4  (10) 

^  -  (1-4A3)/A3  c2  -  (1-4A6)/A6  . 

The  L  x  L  determinant  in  (9)  can  only  be  evaluated  numerically  for  a  given  set 

of  coupling  constants  J^j.  As  all  the  coefficients  A^  involve  through  (6), 

T  can  be  obtained  as  a  function  of  the  J’s.  The  results  of  this  calculation 
c 

are  discussed  in  the  following  section. 

III.  Results  and  discussion 

Throughout  this  paper,  we  take  as  the  unit  of  energy,  and  the  length 
is  measured  in  the  unit  of  the  lattice  constant  in  our  numerical  work.  We 


! 

first  show  the  critical  temperature  as  a  function  of  the  unit  cell  width  L  in  < 
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Fig.  1.  Four  cases  of  different  values  are  plotted  for  fixed  J&  and  Jab. 
We  first  note  that  T  decreases  with  increasing  L  ,  in  general.  This  is 

C  cl 

because  we  have  assumed  J  <  J,  in  our  calculation,  and  it  is  well  known  that 

a  b 

the  bulk  transition  temperature  is  approximately  proportional  to  the  exchange 
coupling  if  the  spin  and  lattice  structure  are  specified.  For  a  given  L  ,  it 

a 

is  observed  that  Tc  increases  steadily  as  L  increases  indefinitely.  All  the 
curves  approach  the  same  limiting  value  Tc  -  5.073  J^/kg,  which  is  simply 

This  is  of  course  easily  understood. 

c 

To  study  the  effects  of  the  interface  exchange  coupling  on  the 

transition  temperature  of  the  superlattice,  T£  is  calculated  for  various  Jab 

but  fixed  J  and  L  .  Figure  3  shows  the  results  for  J  -  J,  .  Curve  c 

corresponds  to  the  case  Jab  -  and  appears  like  that  of  a  slab  of  material 

B.  It  is  very  interesting  to  see  that  there  exists  a  critical  interface 

coupling  J^b  -  1.46  Jb  such  that  Tc  remains  constant  for  any  Lb ,  as 

represented  by  curve  b.  For  J^b  larger  than  this  critical  value,  we  find  that 

T  is  higher  than  both  and  T^  .  Curve  a  illustrates  one  such  particular 

c  c  c 

case  Jflb  -  2Jb.  This  suggests  that  there  exists  an  interface  magnetism  in  the 

system.  For  Jab  >  Jab,  the  system  may  order  in  the  interface  layers  before  it 

orders  in  other  layers.  The  situation  does  not  change  for  J  *  Jg.  The 

results  for  J  -  0.5  J.  are  plotted  in  Fig.  4,  whose  curves  are  very  similar 

except  that  the  critical  interface  coupling  J^b  in  this  case  is  bigger  than  in 

Fig.  3.  This  is  because  the  lover  T^a^  needs  a  stronger,  compensating 

c 

interface  exchange  coupling  in  order  to  reach  the  interface  phase  transition. 

Figure  5  demonstrates  that  as  long  as  -  Jb,  the  critical  coupling  J^b 

-  1.46  Jb  does  not  depend  upon  the  thickness  Lft.  Curve  d  remains  the  same  for 

all  three  cases  calculated,  while  for  J  .  -  2.0  J.  and  0. 1  J.  ,  T  behaves  very 

ab  o  oc 

differently  when  is  varied.  Curves  (a.a*),  <b,b')  and  (c,c')  correspond  to 
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La  -  3,  5  and  10  respectively.  For  J  ^  when  the  interface  magnetism 

appears,  T  decreases  as  L  increases.  For  the  a  and  b  cases,  material  B 
c 

dominates  the  system  for  L  >  13.  For  the  case  c,  ™  10,  and  Tc  drops 
quickly  at  the  beginning  and  approaches  T^^  for  L  >  20.  On  the  other  hand, 
the  transition  temperature  increases  with  L  for  J  ^  <  J^.  The  apparently 
anomalous  behavior  of  curve  c'  can  be  understood  in  the  following  manner.  As 
-  0.1  J^,  Tc  is  mainly  determined  by  material  A(B)  if  1^(1^)  is  larger 
than  L^(La>.  Hence,  changes  little  until  L  approaches  18,  because  L&  -  10 
in  this  case. 

Figure  6  shows  the  phase  diagram  of  the  system  by  plotting  the  critical 
interface  exchange  coupling  as  a  function  of  J  .  Three  different  thicknesses 
are  plotted.  It  is  interesting  to  note  that  the  three  curves  meet  when 
J&  -  J^,  suggesting  that  the  critical  interface  exchange  is  independent  of  the 
relative  thickness  of  the  constituents  in  the  unit  cell  as  long  as  their 
exchange  couplings  are  the  same.  This  is  of  course  in  agreement  with  the 
above  results  in  Fig.  5. 

We  have  considered  ferromagnetic  interface  coupling  in  our  discussion. 
The  critical  temperature  of  the  system,  however,  remains  the  same  if  the 
interface  exchange  coupling  becomes  antiferromagnetic.  This  is  consistent 
with  the  discussion  in  Ref.  7. 

IV.  Conclusions 

We  have  investigated  the  criticality  of  a  two-component  ferromagnetic 
superlattice  by  considering  one  unit  cell,  assuming  nearest- neighbor  spin 
exchange  couplings.  The  interface  exchange  energy  is,  in  general,  different 
from  either  of  the  bulk  ones.  For  the  first  time,  the  critical  value  JaI>  is 
introduced.  When  is  larger  than  this  critical  value,  the  interface 
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3  b 

magnetism  appears,  and  T£  for  the  system  is  higher  than  either  Tc  or  Tc-  When 

3  b 

J  ,  is  smaller  than  the  critical  value,  T  <  T  ,  T  . 
ab  c  c  c 

Finally,  we  remark  that  the  method  is  general  and  can  be  used  to  treat 
the  ferromagnetic/antiferromagnet  and  antiferromagnet/antif erromagnet 
superlattices.  With  modifications  in  spatial  fluctuations  of  the  coupling 
J's,  the  method  can  easily  be  generalized  to  superlattices  involving  certain 
type  of  amorphous  magnetic  materials. 
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Figure  captions 

1.  Unit  cell  of  the  superlattice  composed  of  ferromagnetic  materials  A  and  B, 

where  L  -  is  the  thickness  of  the  cell.  The  solid  and  dashed 

lines  represent  the  nearest-neighbor  exchange  coupling  energies  between 
spins  in  one  material  and  in  the  interface,  respectively. 

2.  Dependence  of  Tc  on  the  thickness  L.  The  interface  coupling  is  chosen  as 
Jak  -  0.9  and  Ja  -  0.8  J^.  Number  of  spin  layers  of  material  A  is  (a) 
3,  (b)  4  and  (c)  5,  (d)  6. 

3.  T  as  a  function  of  L  for  L  -  3  and  J  -  J,  but  various  interface 

c  a  a  b 

coupling  strengths.  (a)  J  ^  -  2.0  J^,  (b)  -  1.46  and  (c  >  Jab  "  V 

4.  Same  as  Fig.  3  except  that  Ja  -  0.5  J^.  (a)  -  3.0  J^,  (b)  -  2.51 

J.  and  (c)  J  ,  -  1.9  J,  . 
b  ab  b 

5.  T  as  a  function  of  L  for  J  -  J.  but  different  J  ,  and  L  .  J  ,  -  2.0  J. 

c  a  b  ab  a  ab  b 

for  curves  a,  b  and  c,  and  J  b  -  0.1  for  a',  b'  and  c' .  The  number  of 
spin  layers  La  is  3  for  a,  a',  5  for  b,  b' ,  and  10  for  c,  c'. 

6.  Phase  diagram  in  terms  of  the  coupling  J  ^  and  J&  for  La  -  (a)  3,  (b)  4 


and  (c)  10. 
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